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Abstract 

We express a matrix version of the self-induced transparency (SIT) equations in the bidiffer- 
ential calculus framework. An infinite family of exact solutions is then obtained by application 
of a general result that generates exact solutions from solutions of a linear system of arbitrary 
matrix size. A side result is a solution formula for the sine-Gordon equation. 

1. Introduction. The bidifferential calculus approach (see [I] and the references therein) aims to 
extract the essence of integrability aspects of integrable partial differential or difference equations 
(PDDEs) and to express them, and relations between them, in a universal way, i.e. resolved from 
specific examples. A powerful, though simple to prove, result [l]-[3] (see section 6) generates families 
of exact solutions from a matrix linear system. In the following we briefly recall the basic framework 
and then apply the latter result to a matrix generalization of the SIT equations. 

2. Bidifferential calculus. A graded algebra is an associative algebra f2 over C with a direct 
sum decomposition Q, = (J) r>0 Q r into a subalgebra A := 0° and „4-bimodules f2 r , such that 
ST £l s C rt r+s . A bidifferential calculus (or bidifferential graded algebra) is a unital graded algebra 
f2 equipped with two (C- linear) graded derivations d, d : f2 — > ft of degree one (hence dSl r C $7 r+1 , 
df2 r C f2 r+1 ), with the properties 

d 2 z = VzeC, where d z := d - z d , (1) 

and the graded Leibniz rule d z (xx') = (d 2 x) x' + ( — l) r X&zX\ f° r ai l X G ^ r an d x' & 

3. Dressing a bidifferential calculus. Let (£1, d, d) be a bidifferential calculus. Replacing d z 
in ([!]) by D z := d — A — zd with a 1-form A £ ft 1 (in the expression for D z to be regarded as a 
multiplication operator), the resulting condition = (for all z G C) can be expressed as 

dA = = dA- AA . (2) 

If ^ is equivalent to a PDDE, we have a bidifferential calculus formulation for it. This requires 
that A depends on independent variables and the derivations d, d involve differential or difference 
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operators. Several ways exist to reduce the two equations ^ to a single one: 

(1) We can solve the first of ^ by setting A = dcf>. This converts the second of ^ into 

dd<f> = dcf) d<f> . (3) 

(2) The second of ^ can be solved by setting A = (dg) g -1 . The first equation then reads 

d({dg)g- 1 )=Q. (4) 



(3) More generally, setting A = [dg — (dg)A\ g^ 1 , with some A E A, we have dA — A A = 
(dA) gAg^ 1 + (dg) (dA — (dA) A) g^ 1 . As a consequence, if A is chosen such that dA = (dA) A, 
then the two equations ^ reduce to 



d([dg-(dg)A]g-') = 0. (5) 

With the choice of a suitable bidifferential calculus, ^ and Q, or more generally ([5]), have been 
shown to reproduce quite a number of integrable PDDEs. This includes the self-dual Yang-Mills 
equation, in which case ([3]) and Q correspond to well-known potential forms [Tj. Having found a 
bidifferential calculus in terms of which e.g. ([3]) is equivalent to a certain PDDE, it is not in general 
guaranteed that also Q represents a decent PDDE. Then the generalization §5§ has a chance to 
work (cf. |lj). In such a case, the Miura transformation 

[dg - (dg)A] g- 1 = d<p (6) 

is a hetero-Backlund transformation relating solutions of the two PDDEs. 

Backhand, Darboux and binary Darboux transformations can be understood in this general 
framework [I], and there is a construction of an infinite set of (generalized) conservation laws. 
Exchanging d and d leads to what is known in the literature as 'negative flows' [jjj. 

4. A matrix generalization of SIT equations and its Miura-dual. A = Mat(n, n, C°°(M 2 )) 
denotes the algebra of n x n matrices of smooth functions on M 2 . Let Q. = A 8) A(C 2 ) with the 
exterior algebra /\(C 2 ) of C 2 . In terms of coordinates x,y of M 2 , a basis Cl>C2 of /\ 1 (C 2 ), and a 
constant n x n matrix J, maps d and d are defined as follows on A, 

df= 1 2[JJ]®(l + fy®(2, df = f x ®tl + \[J,f]®C2 

(see also El). They extend in an obvious way (with dQ = d(i = 0) to Cl such that (fi, d, d) becomes 
a bidifferential calculus. We find that ^ is equivalent to 

^ = |[[J,^],^-ij] . (7) 

Let n = 2m and J = block-diag(I, — I), where I = I m denotes the m x m identity matrix. 
Decomposing <p into m x m blocks, and constraining it as follows, 

P q ), (8) 
Q -PJ 

([7]) splits into the two equations 

Pxy = (q 2 ) y , q xy = q- p y q - qp y ■ (9) 
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We refer to them as matrix-SIT equations (see section 5), not purporting that they have a similar 
physical relevance as in the scalar case. The Miura transformation Q (with A = 0) now reads 

-i 



9x9 



2 [ J >' 



W,9]g- 1 



'y ■ 



(10) 



Writing 



a b 
c d 



with m x m matrices o, b, c, d, and assuming that a and its Schur complement S(a) = d — ca 1 b is 
invertible (which implies that g is invertible), (10) with ^ requires 

b = —ca~ 1 d, a x = —c x a~ 1 c, d x = —c x aT l ca d . (11) 

The last equation can be replaced by d x d~ l = a^a -1 . Invertibility of S(a) implies that d and 
I + r 2 are invertible, where r := ca -1 . The conditions (11) are necessary in order that the Miura 
transformation relates solutions of ^ to solutions of its 'dual' 

{9x9^% = \ \gJg~ x ,J\ , 

obtained from Q. Taking (11) into account, the Miura transformation reads 

q = -c x a~ x = -r x -ra x a~ x , q y = —r (I + r 2 ) -1 , p y =/-(/ + r 2 ) -1 . 

As a consequence, we have 
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% +Py 



Py ■ 



(12) 
(13) 
(14) 



Furthermore, the second of ( 11 ) and the first of ( 13 ) imply a x a 1 = qr. Hence we obtain the system 

-r(/ + r 2 )-\ 



r qr . 



Qy 



(15) 



which may be regarded as a matrix or 'noncommutative' generalization of the sine-Gordon equation. 
There are various such generalizations in the literature. The first equation has the solution q = 
— X^feLo( — rk Tx r * ' i ^ the sum exists. Alternatively, we can express this as q = — (I+rLrn)~ l (r x ), 
where (tr) denotes the map of left (right) multiplication by r. This can be used to eliminate q 
from the second equation, resulting in 



{{I + r L r R )-\r x )) = r {I + r 2 )' 



(16) 



If r = tan(0/2) 7r with a constant projection 7r (i.e. 7r 2 = tt) and a function 6, then (16) reduces 
to the sine- Gordon equation 



'XL/ 



sin ( 



(17) 



(| 1 5 1) can be obtained directly from \12\ as follows, by setting 

-l 
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I 



hence g 



-i 



I r 
-r I 



(I + r 2 



This leads to 

((r x r + rpr + p)(I + r 2 )^ 1 )^ 



0, ((r x +rp - pr)(I + r 2 ) l ) t = r(I + r 2 )" 1 



where p := a x a l . Setting an integration 'constant' to zero, the first equation integrates to p 



r x r 



(ra x a 



rpr. With its help, the second can be written as {r x + rp) y = r(I + 



Since q 



r p, this is the second of (15). The first follows noting that qr = p. 
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5. Sharp line SIT equations and sine-Gordon. We consider the scalar case, i.e. m = 1. 
Introducing £ = 2^[aq with a positive constant a, V = 2q y , M = 2p y — 1, and new coordinates z, t 
via x = ^/a{z — t) and y = y/az, the system Q is transformed into 



V t = £M, M t = -£V, 
and the relation between £ and V takes the form 

£ z + £ t = aV . 

These are the sharp line self-induced transparency (SIT) equations (5|j7j. We note that V 2 + M 2 
is conserved. Indeed, as a consequence of (14), we have V 2 + N 2 = 1. Writing V = — sm9 and 



M = — cos 9, reduces the first two equations to £ = 8t- Expressed in the coordinates x, y, the third 



then becomes the sine-Gordon equation (17) (cf. 16]). As a consequence of the above relations, q 
and p depend as follows on 6, 

q = -\0x, q y = -\sin6, p y = |(1 - cos#) . (18) 



These are precisely the equations that result from the Miura transformation (10) (or from (13)), 
choosing 
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COS 2 




and (12) becomes the sine-Gordon equation (17). The conditions (11) are identically satisfied as a 



consequence of the form of g. 

6. A universal method of generating solutions from a matrix linear system. 

Theorem 1. Let (f2,d, d) be a bidifferential calculus with O = A.® /\(C 2 ), where A is the algebra 
of matrices with entries in some algebra B (where the product of two matrices is defined to be 
zero if the sizes of the two matrices do not match). For fixed N,N', let X G Mat(iV, N, B) and 
Y G Mat(iV, N, B) be solutions of the linear equations 

dX = (dX)P, dY = (dY)P, RX - X P = -QY , 

with d- and d-constant matrices P,R G Mat(iV, N, B), and Q = V U , where U G Mat(n, N',B) 
and V G M&t(N,n, B) are d- and d-constant. If X is invertible, the nx n matrix variable 

cj) = UYX~ X V G Mat(n,n,B) 

solves d<p = (dcj>) <f> + di9 with = UY X 1 RV , hence (by application of d) also Q). □ 

There is a similar result for ^ |3j. The Miura transformation is a corresponding bridge. 

7. Solutions of the matrix SIT equations. From Theorem [T] we can deduce the following 
result, using straightforward calculations |8|, analogous to those in [2] (see also [3]). 

Proposition 2. Let S G Mat(M, M, C) be invertible, U G Mat(m,M,C) ; V G Mat(M,m,C), and 
K G Mat(M, M, C) a solution of the Sylvester equation 

SK + KS = VU . (19) 

Then, with H = e~ Sx ~ s 1 v and any po G Mat(m, m, C) (more generally x - dependent) , 

q = US(I M + (KS) 2 )- 1 V, p = Po -USKS(I M + (KS) 2 )- 1 V (20) 

(assuming the inverse exists) is a solution of |9p. □ 



4 



If the matrix S satisfies the spectrum condition 



a(S)na(-S) 



(21) 



(where o~(S) denotes the set of eigenvalues of S), then the Sylvester equation (19) has a unique 

3.g. [9]. 



see e. 



solution K (for any choice of the matrices U, V) 

By a lengthy calculation [8] one can verify directly that the solutions in Proposition [2] satisfy 
Alternatively, one can show that these solutions actually determine solutions of the Miura 



(14) 



transformation (cf. [3]), and we have seen that (14) is a consequence. 

There is a certain redundancy in the matrix data that determine the solutions (20) of Q. This 
can be narrowed down by observing that the following transformations leave (19) and (20) invariant 
(see also the NLS case treated in [2]). 

(1) Similarity transformation with an invertible M G Mat(M, M, C) : 



MSM 



-i 



K ^ MKM 



-i 



V ^ MV , U^UM 



-i 



As a consequence, we can choose S in Jordan normal form without restriction of generality. 
(2) Reparametrization transformation with invertible A,B G Mat(M, M, C) : 



K ^ B X KA 



V 



B~ X V 



U i->- UA 



-i 



S i — y ABB 



(3) Reflexion symmetry: 

St-t-S, K^-K 1 , V^K^V, U^UK- 1 , p i-> Po - UK X V . 

This requires that K is invertible. More generally, such a reflexion can be applied to any Jordan 
block of S and then changes the sign of its eigenvalue [8] (see also [2||10| ) . The Jordan normal form 
can be restored afterwards via a similarity transformation. 
The following result is easily verified [8j. 

Proposition 3. Let S, U, V be as in Propositions and T G Mat(M, M, C) invertible. 



(1) Let T be Hermitian (i.e. =T) and such that 



TST 

-l 



U = V ] T. Let K be a solution 



of (19), which can then be chosen such that = TKT L . Then q and p given by (20) with 
Po are both Hermitian and thus solve the Hermitian reduction of 



(2) Let T = T~ L (where the bar means complex conjugation) and S = TST 1 , U = UT 1 and 
V = TV. Let K be a solution of (19), which can then be chosen such that K = TKT 1 . Then q 
and p given by (20) with po = po satisfy q = q and p = p, and thus solve the complex conjugation 
reduction of |$|). □ 

8. Rank one solutions. Let M = 1. We write S = s, U = u, V = v J , K = k (where T means 
the transpose) and S = £ = e" 



sx — s 



. Then (19) yields k = (v J u)/(2s). From (20) we obtain 



2sk£ 
1 + (k0 [ 



IT, P=P0 + 



2s 



i + (ko- 



7T, po := po ~ 2s TV , TV 



The Miura transformation (13) implies r = — q y (I — p y ) , and we obtain 

2fc£ 

r --l-{ki)^' 
which is sing ular. But 9 = -2 arctan(2/c£/[l - (fe£) 2 ]) is the 

single kink solution of the sine-Gordon 



equation (17). 
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9. Solutions of the scalar (sharp line) SIT equations. We rewrite p in (20), where now 
m = 1, as follows, 



p = po - tr ((SK + KS) SKS (I M + (KH) 2 ) -1 ) 
= p + tr ((/ M + (KS) 2 ), (I M + (^H) 2 )" 1 ) 
= p + (logdet (Jm + ^S) 2 ))^ , 



(22) 



using (19) and the identity (det M) x = tr(M x M 1 ) detiW for an invertible matrix function M. 
q in ( 20 ) can be expressed as 



q = 2ti(SKS(I M + (KS) 2 )- 1 ) . 
In particular, if S is diagonal with eigenvalues Sj, i = 1, . . . , M, and satisfies (|21|), then the solution 



K of the Sylvester equation (19), which now amounts to r&nk(SK + KS) = 1, is the Cauchy-type 
matrix with components Kij = v\ Uj/ (s{ + sj), where u%,Vi 6 C. Figs. [T] and [2] show plots of two 



examples from the above family of solutions. 





Figure 1: A scalar 2-soliton solution with S = diag(l, 2) and Ui = Vi = 1. 




10. A family of solutions of the real sine-Gordon equation. Via the Miura transformation 
(18), Proposition [2] determines a family of sine-Gordon solutions (see also e.g. [6,11-16 for related 
results obtained by different methods). 

Proposition 4. Let S G Mat(M, M, C) be invertible and K £ Mat(M, M, C) such that r&nk(SK+ 
KS) = 1, det(I M + {Ka) 2 ) G R with 3 = e " Sa; - s ~\ and tv{SKa (J M + (KS) 2 )- 1 ) iM (Wiere 
i is i/ie imaginary unit). Then 



fl = 4arctan( i + ^_^ ) urift /3 := (log | det(/ M + (KH) 2 )))^ (23) 
solves the sine- Gordon equation 9 xy = sin 9 in any open set ofR 2 where det(I M + (KS) 2 ) ^ 0. 
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Proof: Let p be given by (22). Due to the assumption det(/M + {Ka) 2 ) G M, p y is real, hence 
(14) implies |1 — 2p y \ 2 = 1 — ^q y 2 . It follows that q y 2 is real. Since another of our assumptions 
excludes that q y is imaginary, it follows that |1 — 2p y \ < 1. Hence the equation cos 6 = 1 — 2p y 
(second of (18)) has a real solution 9. Inserting the expression (22) for p, we arrive at cos 9 = 
1-2 (logdet(iM + (KS) 2 )) x . Moreover, (14) shows that p y > and thus < p y < 1. Using 
identities for the inverse trigonometric functions, we find (23), where f3 = p y . □ 

Proposition [3] yields sufficient conditions on the matrix data for which the last two assumptions 
in Proposition [4] are satisfied. 
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